Mobility in semiconducting carbon nanotubes at finite carrier density 
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Carbon nanotube field-effect transistors operate over a wide range of electron or hole density, 
controlled by the gate voltage. Here we calculate the mobility in semiconducting nanotubes as a 
function of carrier density and electric field, for different tube diameters and temperature. The 
low-field mobility is a non-monotonic function of carrier density, and varies by as much as a factor 
of 4 at room temperature. At low density, with increasing field the drift velocity reaches a maximum 
and then exhibits negative differential mobility, due to the non-parabolicity of the bandstructure. 
At a critical density p c ~ 0.35-0.5 electrons/nm, the drift velocity saturates at around one third of 
the Fermi velocity. Above p c , the velocity increases with field strength with no apparent saturation. 



Carbon nanotubes have enabled new nano-electronic 
devices such as quasi-one-dimensional field-effect tran- 
sistors (CN-FET) i] and electro-optical devices H 0. 
Ballistic devices have generated excitement as the ideal 
limit of transistor performance Yet diffusive devices 
remain important 0, H °-ue to the very high mobil- 
ity even at room temperature 0. Much effort has gone 
into determining their transport properties experimen- 
tally HH USED, 03 and theoretically [H Q El 03 

Conventional transistors operate by control of carrier 
density, i.e. band filling via an external gate electrode 
and its capacitive coupling to the conduction channel. 
Yet while carrier density is a crucial variable in device 
operation, relatively little is known about how the mo- 
bility in nanotube transistors varies with carrier density, 
especially at high fields. Only recently have systematic 
studies of the dependence on gate voltage been reported 

mm 

Here we calculate the mobility in carbon nanotubes as 
a function of charge carrier density. We find that the mo- 
bility depends sensitively on carrier density. Specifically, 
at low field and low temperature, the mobility initially 
increases with carrier density, because the phase space 
for scattering mi is reduced by band filling. The mobil- 
ity then decreases abruptly as the band filling reaches the 
second band, because a new channel for scattering opens 
up. The low-field mobility can be expressed as the prod- 
uct of scattering time and hinverse effective mass, where 
the effective mass is nearly independent of temperature, 
and increases smoothly with band filling. Our results 
quantify a recent analysis of low- field mobility . 

At high fields, we find negative differential mobility 
(NDM) at low density, decreasing and disappearing with 
increasing density. Surprisingly, this NDM appears to 
be unrelated to the occupation of the second band, but 
can be ascribed to the nonparabolicity of the first band. 
The minimum differential mobility, i.e. maximum mag- 
nitude of the NDM, increases with diameter and de- 
creases with temperature, with a systematically larger 
value in tubes with mod(n — m, 3) = 1 than in those 
with mod(n — to, 3) = — 1. This is due to the stronger 
band nonparabolicity in the former. 



Our calculations use a tight-binding description for 
the electronic structure The phonons are described 
using a model similar to that of Aizawa et al. 
We model the electron-phonon interaction by the Su- 
Schrieffer-Heeger model |2(1 ]. with matrix element t = 
to — gSu dependent on the change of the nearest neighbor 
C-C distance (Su), where t = 3 eV. We take the electron- 
phonon coupling constant to be g — 5.3 eV/A |2lj| . 

We solve the steady-state multi-band Boltzmann equa- 
tion in the presence of an electric field, to determine the 
electron distribution function g k : 
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where Wk,k+q is an electron-phonon scattering rate |l3| . 
The indices k and q label both the continuous wave vec- 
tor along the tube axis and the band index. Using the 
non-equilibrium distribution function g k from Eq. , we 
calculate the drift velocity Vd and the mobility fi accord- 
ing to: 
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where V k is the band velocity and gt = ^2i k 9k 'is the band 
filling factor (electrons/carbon). 

The low-field electron mobility is shown in Fig. for 
a (19,0) tube. (The hole mobility is identical, due to 
the electron-hole symmetry in our model.) The mobil- 
ity has a complex and highly non-monotonic dependence 
on density at low temperature, and decreases at higher 
temperature due to phonon scattering. 

To understand the low-field mobility, we use the clas- 
sical relationship [i = er /m av to decompose the mobility 
into a product of two terms: the scattering time t and 
the averaged inverse effective mass m~ t ,. The natural 
definition of m av and resulting definition of r are |22| : 
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FIG. 1: (a) Low-field electron mobility in (19,0) tube (d=1.5 
nm) as a function of the charge carrier density p. (b-c) De- 
composition of mobility into a product of (b) inverse effective 
mass and (c) scattering time. Temperatures are (from top to 
bottom in each frame) T=150 K (dashed green), T=200 K 
(dashed blue), T=250 K (dashed red), T=300 K (solid ma- 
genta), T=350 K (solid green), T=400 K (solid blue), T=450 
K (solid red). The vertical dashed line is at p = pd — 0.99 
electrons/nm (see text). The density of 1 electrons/nm corre- 
spond to band filling of 0.0057 electrons/carbon for a (19,0) 
tube. 



The effective mass contribution to the low-field mobility 
is shown in Fig. \\]p - it is virtually independent of tem- 
perature and it is a monotonically decreasing function of 
the electron density p. In contrast, the scattering time 
contribution in Fig. ^ has a strong dependence on both 
temperature and charge density. The charge scale is de- 
termined primarily by the density pd at which the first 
band is filled up to the energy of the bottom of the second 
band (at T = K). The scattering time drops sharply 
around pd, reaching a minimum before rising again. 

At low field, g k is given essentially by the Fermi-Dirac 
distribution, and we can solve Eq. for m av analytically 
for a given band dispersion . We use a hyperbolic band 
dispersion 
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which provides a good approximation to the tight- 
binding bands in semiconducting carbon nanotubcs. 
Here A is half the bandgap energy and too is the effective 
mass at the bottom of the band. At zero temperature, 
we solve Eq. (j3J and Eq. (@J for the density-dependent 
effective mass: 
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where the relationship Kp — up/ A between the ID Fermi 
vector Kp and density p has been used. As the Fermi level 
increases with density, the effective mass gradually in- 
creases and almost doubles when the Fermi level reaches 
the bottom of the second band (see Fig.^Jj). 

With increasing density, the scattering time initially 
increases. This can be understood as resulting from the 
decreasing density of states available for scattering. With 
further density increase, electrons start to populate the 
second band and an additional channel opens up for scat- 
tering. This reduces the scattering time at all tempera- 
tures, with the largest change occurring at low temper- 
ature. As the temperature increases, the density depen- 
dence of the scattering time becomes weaker. Therefore, 
above room temperature, both the effective mass and 
the scattering time make comparable contributions to the 
density dependence of the low-field mobility. In contrast, 
at low temperatures, the scattering time dominates the 
density dependence of the low-field mobility. 

The non-monotonic density dependence of the mobil- 
ity, predicted here, agrees well with experiment 
Experiment reports a maximum in the mobility and a 
few-fold drop from its peak value as the gate voltage 
increases by several volts. This behavior is nicely re- 
produced in Fig. n Agreement in the absolute value of 
the mobility is not expected, because we only consider 
phonon scattering, while additional scattering mecha- 
nisms (e.g. impurities, charged defects in the oxide) are 
present in the experiment. 

Carbon nanotube devices can be made with a range of 
nanotube diameters, offering different bandgaps. We find 
the scattering time to be approximately proportional to 
the diameter for all charge densities and temperatures. 
Indeed, the combination r/d depends primarily on tem- 
perature and charge density and only weakly on diameter 
as shown in Fig. [2 (In the zero density limit, there is an 
additional simplification 

So far we have focused on the transport properties at 
low fields. Transport at high fields is also important, 
and our results for the dependence of carrier velocity on 
electric field are summarized in Fig. For low carrier 
density, the drift velocity has a maximum at a field of 
E ~ 0.5-1.0 (V/jLtm), which is within the range of fields 
used in normal device operation, and a negative differen- 
tial mobility (NDM) at higher fields. The NDM becomes 
weaker with increasing density, and disappears entirely 
at a critical density p c . We find that p c is not a strong 
function of tube diameter d; for the range of d considered 
here, p c ~0.35-0.5 electrons/nm at room temperature. 
Above pc, the drift velocity Vd is a monotonic function of 
electric field, with no saturation behavior up to the fields 
studied here |23j . 

We can fit all the results in Fig. [3^ using the following 
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FIG. 2: Diameter scaling of the scattering time in carbon 
nanotubes. We show scaled scattering time r/d versus scaled 
density p/ pd for zig-zag (n,0) tubes with diameter range d = 
1.0 - 2.0 nm. In electrons/nm, pd = 1.44 for n=13 (solid 
red), pd = 1-33 for 14 (solid blue), pd = 1.17 for 16 (solid 
green), pd = 1.10 for 17 (solid magenta), pd ~ 0.99 for 19 
(solid black), p d = 0.94 for 20 (dashed red), p d = 0.85 for 
22 (dashed blue), pd = 0.81 for 23 (dashed green), pd = 0.75 
for 25 (dashed magenta) at temperatures (a) T = 150 K, (b) 
T = 300 K, (c) T = 450 K. 



phenomenological formula: 



V d = 
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where p,Q is the low-field mobility (from Fig.^), and V s2 is 
the saturation velocity at p c . The value of V s2 = 3 x 10 7 
(cm/s) is about 1/3 of the Fermi velocity for graphene 
and it is roughly independent of tube diameter. The 
other two parameters, V s ± and E 2 , are fitted separately 
for each density. 

A recent experiment 0] a ^ high drain voltage in CN- 
FETs (with d ~ 2.0-2 .4 nm) gave a saturation velocity of 
2 x 10 7 cm/s. This is a factor of 2.5 smaller than the max- 
imum velocity predicted in the zero density limit fl4ll2l| . 
The electrostatic gate capacitance per unit length was es- 
timated to be C g = 2 x 10~ n F/m. This gives a charge 
density p = C g V g = 1.1 electrons/nm for the gate voltage 
of 9 V [12J . This implies that the experimental device was 
in the regime, where the carrier density p was well above 
critical density p c . Therefore, the zero density limit is 
not applicable, and the finite density result obtained from 
Fig. agrees much better with the experiment. 

Negative differential mobility has been predicted theo- 
retically 0,0,0] and negative differential conductance 
has been observed experimentally |25l|. One possible in- 
terpretation of this behavior [I3. ll5| is that the occupa- 
tion of the second band (which has higher effective mass) 



FIG. 3: (a) Drift velocity (filled circles) at T = 300 K, 
as a function of electric field, in (19,0) tube (d =1.5 nm) 
at different charge densities. Curves show best fit to Eq. ED 
Charge densities (in electrons/nm) are, from top to bottom, 
p = limit (red), p = 0.25 (blue), 0.5 (green), 0.75 (red), 
1.0 (blue), 1.5 (green), 2.0 (blue), (b) The p = results 
from (a) are repeated here (red solid curve, filled circles), 
and compared with different simplified models: single band 
(black open squares); single hyperbolic band (black dotted 
curve); and single parabolic band (red open circles), extrap- 
olated as red dashed curve using the best fit to Eq. JSJ , i.e. 
V a i = 5.0 10 7 cm/s and E2 = 00. 



is responsible for the NDM. In agreement with previous 
studies 0, 0| , we find that at high field a significant 
fraction of charge carriers occupy the second band. How- 
ever, we find that this alone does not explain the NDM. 
In fact, if we artificially restrict the electrons to the first 
band, the NDM becomes stronger rather than weaker 
(Fig. Eh). 

To understand the origin of the NDM, we examine the 
effect of changes in bandstructure. The tight-binding 
bandstructure is well described by a hyperbolic disper- 
sion [Eq. Q]; and Fig. [3Jd shows that the hyperbolic 
bandstructure gives results essentially identical to the 
real bandstructure (when only the first band is included) . 
However, when we replace the realistic hyperbolic band- 
structure with a parabolic band in Fig. EJd, the NDM 
disappears entirely. Therefore, we conclude that the non- 
parabolicity of the first band is responsible for the NDM 
in semiconducting carbon nanotubes. 

In the NDM regime, the mobility is approximately con- 
stant over a wide range of electric field, and equal to the 
most negative value, i.e. the maximum NDM. Therefore, 
we can characterize the NDM of a given tube by two pa- 
rameters: the most-negative value p, n ; and the field at 
which the velocity reaches a maximum. 

With these parameters, we can summarize the behav- 
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FIG. 4: (a) Maximum NDM, function of charge 

density in (19,0) tube for different temperatures from T = 450 
K (top) to T = 150 K (bottom) with in steps of 50 K. (b) \x n 
as a function of charge density in zig-zag tubes at T = 300 
K. The mod(n, 3) = 1 tubes are shown by the solid curves 
n = 16 (blue), 19 (green), 22 (magenta), 25 (black), and 
mod(n, 3) = —1 tubes are by the dashed curves n=17 (blue), 
20 (green), 23 (magenta), (c) The electric field corresponding 
to the maximum drift velocity in the same tubes as in (b). 



ior over a range of density, different temperatures and 
tube diameters as shown in Fig.0] As the temperature is 
increased, the maximum NDM decreases approximately 
linearly with temperature. 

The maximum NDM, shown in Fig.^D, increases with 
tube diameter and has a clear family effect. The family 
of tubes with mod(n — m, 3) = 1 have higher NDM than 
family of tubes with mod(n — to, 3) = —1. This can be 
understood by the fact that the ratio A/too measures 
the non-parabolicity of the band, which depends on the 
tube family. In the calculations for a single hyperbolic 
band Eq. |0J, we find that the NDM is a linear function of 
A/too, even though A and mo are independent variables. 
If A and too coincide with the tight-binding values, the 
drift velocity is virtually the same as in the tight-binding 
band calculations restricted to a single band (as show in 

Fig. Eh). 

Finally, the velocity reaches a maximum at a field E m 
shown in Fig.QJ;, which decreases with the tube diameter. 
This suggests that conditions for the NDM are easier to 
achieve in larger diameter tubes. At the critical density 
p c , the field E m grows to infinity and NDM disappears. 

In conclusion, we calculated the transport properties 
of CN-FETs as a function of change density as controlled 



by the gate potential. The low-field mobility dependence 
on the charge density is mainly determined by the scat- 
tering at low temperatures, and by both scattering and 
non-parabolicity of the first conduction band at high tem- 
peratures. At high fields, the drift velocity shows a nega- 
tive differential mobility behavior which is large in large 
diameter tubes and low temperatures. We showed that 
the non-parabolicity of the first conduction band is re- 
sponsible for the NDM behavior. 
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